Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis     

Ramin Pakzad  Shanyong Wang  Scott W. Sloan 《国际地质力学数值与分析法杂志》2020,44(10):1472-1494

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.  相似文献   

Effects of the geometry of two‐dimensional fractures on their hydraulic aperture and on the validity of the local cubic law          下载免费PDF全文

Ahmad Hajjar  Luc Scholtès  Constantin Oltéan  Michel Antoine Buès 《水文研究》2018,32(16):2510-2525

Flow through rough fractures is investigated numerically in order to assess the validity of the local cubic law for different fracture geometries. Two‐dimensional channels with sinusoidal walls having different geometrical properties defined by the aperture, the amplitude, and the wavelength of the walls' corrugations, the corrugations asymmetry, and the phase shift between the two walls are considered to represent different fracture geometries. First, it is analytically shown that the hydraulic aperture clearly deviates from the mean aperture when the walls' roughness, the phase shift, and/or the asymmetry between the fracture walls are relatively high. The continuity and the Navier–Stokes equations are then solved by means of the finite element method and the numerical solutions compared to the theoretical predictions of the local cubic law. Reynolds numbers ranging from 0.066 to 66.66 are investigated so as to focus more particularly on the effect of flow inertial effects on the validity of the local cubic law. For low Reynolds number, typically less than 15, the local cubic law properly describes the fracture flow, especially when the fracture walls have small corrugation amplitudes. For Reynolds numbers higher than 15, the local cubic law is valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitudes, and a moderate phase lag between its walls.  相似文献   

Analytical study of ground motion caused by seismic wave propagation across faulted rock masses          下载免费PDF全文

J.C. Li  N.N. Li  S.B. Chai  H.B. Li 《国际地质力学数值与分析法杂志》2018,42(1):95-109

Studying seismic wave propagation across rock masses and the induced ground motion is an important topic, which receives considerable attention in design and construction of underground cavern/tunnel constructions and mining activities. The current study investigates wave propagation across a rock mass with one fault and the induced ground motion using a recursive approach. The rocks beside the fault are assumed as viscoelastic media with seismic quality factors, Q_p and Q_s. Two kinds of interactions between stress waves and a discontinuity and between stress waves and a free surface are analyzed, respectively. As the result of the wave superposition, the mathematical expressions for induced ground vibration are deduced. The proposed approach is then compared with the existing analysis for special cases. Finally, parametric studies are carried out, which includes the influences of fault stiffness, incident angle, and frequency of incident waves on the peak particle velocities of the ground motions.  相似文献   

The bivalve Lembulus bicuspidatus may enhance denitrification in shelf sediment at the Angola–Benguela Frontal Zone     

A Neumann  A Flohr 《African Journal of Marine Science》2018,40(1):91-96

We collected living individuals of the bivalve Lembulus bicuspidatus, which shows an unusual preference for the oxygen-deficient habitat found at the Angola–Benguela Frontal Zone of the southeastern Atlantic. With a series of incubation experiments with ¹⁵N-labelled nitrate as a tracer in combination with membrane-inlet mass spectrometry, we studied the potential contribution of L. bicuspidatus to nitrate reduction in the upper sediment layer. Our preliminary results suggest that L. bicuspidatus enhances nitrate reduction if the oxygen concentration is sufficiently low. The Lembulus-mediated nitrate reduction rate is then similar to the rate of microbial nitrate reduction in the surrounding sediment.  相似文献   

1970-2016年青藏高原岗扎日冰川变化与物质平衡遥感监测研究     

张震  刘时银 《地球信息科学学报》2018,20(9):1338-1349

可可西里处于青藏高原腹地,是青藏高原自然环境的交接与过渡地带。近年来该区域冰川物质平衡可能有从西向东由正转负的趋势,但是其过渡地带岗扎日地区冰川状态未知。本研究利用地形图、SRTM、ASTER和Landsat等资料分析了岗扎日地区冰川面积变化和物质平衡变化,并对可可西里地区冰川变化空间规律进行了探讨,结果表明：①1970-2016年岗扎日冰川总面积年均缩小率为0.08±0.02%。2006年后冰川退缩趋势减缓。②1970-2012年岗扎日冰川平均减薄-8.64±0.30 m,体积减少1.45±0.06 km³,平均物质平衡为-0.21±0.01 m w.e. a^-1。冰川物质平衡趋势由负转正（1970-1999年：-0.34±0.01 m w.e. a^-1;1999-2012：0.16±0.02 w.e. a^-1）。③东南、南、西南朝向作为迎风坡,1970年以来其冰川物质亏损较小,1999-2012年呈现强烈的正平衡。冰川面积变化滞后于物质平衡变化,东朝向和东南朝向冰川面积缩小率最大,主要是因为冰川冰舌较长,末端所处的海拔较低。④气温升高是岗扎日冰川1970-1999年呈现负物质平衡状态的主因,降水增多是1999-2012年正平衡状态的主因。⑤可可西里地区冰川1970s以来面积年均缩小率从西向东不断增大、物质平衡下降,与西风环流和季风环流相关,但局地气候也影响冰川变化和物质平衡。  相似文献   

Explicit stress–strain equations for modeling frictional materials     

Kenneth J. Xu  Martin D. Liu  Buddhima Indraratna 《Marine Georesources & Geotechnology》2018,36(6):722-734

In this paper, a comprehensive study on simulating the shearing behavior of frictional materials is performed. A set of two explicit equations, describing the relationship among the shear stress ratio and the distortional strain and the volumetric strain, are formulated independently. The equations contain three stress parameters and three strain parameters and another parameter representing the nonuniformity of stress and strain during softening. All the parameters have clear physical significance and can be determined experimentally. It is demonstrated that the proposed equations have the capacity of simulating the complicated shearing behavior of many types of frictional materials including geomaterials. The proposed equations are used to simulate the stress–strain behavior for 27 frictional materials with 98 tests. These materials include soft and stiff clays in both reconstituted and structured states, silicon sands and calcareous sands, silts, compacted fill materials, volcanic soils, decomposed granite soils, cemented soils (both artificially and naturally cemented), partially saturated soils, ballast, rocks, reinforced soils, tire chips, sugar, wheat, and rapeseed. It has been demonstrated that the proposed explicit constitutive equations have the capacity to capture accurately the shearing behavior of frictional materials both qualitatively and quantitatively. A study on model parameters has been performed.  相似文献   

Solution for a plane strain rough‐walled hydraulic fracture driven by turbulent fluid through impermeable rock          下载免费PDF全文

N. Zolfaghari  E. Dontsov  A. P. Bunger 《国际地质力学数值与分析法杂志》2018,42(4):587-617

The impact of turbulent flow on plane strain fluid‐driven crack propagation is an important but still poorly understood consideration in hydraulic fracture modeling. The changes that hydraulic fracturing has experienced over the past decade, especially in the area of fracturing fluids, have played a major role in the transition of the typical fluid regime from laminar to turbulent flow. Motivated by the increasing preponderance of high‐rate, water‐driven hydraulic fractures with high Reynolds number, we present a semianalytical solution for the propagation of a plane strain hydraulic fracture driven by a turbulent fluid in an impermeable formation. The formulation uses a power law relationship between the Darcy‐Weisbach friction factor and the scale of the fracture roughness, where one specific manifestation of this generalized friction factor is the classical Gauckler‐Manning‐Strickler approximation for turbulent flow in a rough‐walled channel. Conservation of mass, elasticity, and crack propagation are also solved simultaneously. We obtain a semianalytical solution using an orthogonal polynomial series. An approximate closed‐form solution is enabled by a choice of orthogonal polynomials embedding the near‐tip asymptotic behavior and thus giving very rapid convergence; a precise solution is obtained with 2 terms of the series. By comparison with numerical simulations, we show that the transition region between the laminar and turbulent regimes can be relatively small so that full solutions can often be well approximated by either a fully laminar or fully turbulent solution.  相似文献   

Micromechanical approach to effective viscoelastic properties of micro‐fractured geomaterials          下载免费PDF全文

Cássio B. Aguiar  Samir Maghous 《国际地质力学数值与分析法杂志》2018,42(16):2018-2046

The aim of this paper is to formulate a micromechanics‐based approach to non‐aging viscoelastic behavior of materials with randomly distributed micro‐fractures. Unlike cracks, fractures are discontinuities that are able to transfer stresses and can therefore be regarded from a mechanical viewpoint as interfaces endowed with a specific behavior under normal and shear loading. Making use of the elastic‐viscoelastic correspondence principle together with a Mori‐Tanka homogenization scheme, the effective viscoelastic behavior is assessed from properties of the material constituents and damage parameters related to density and size of fractures. It is notably shown that the homogenized behavior thus formulated can be described in most cases by means of a generalized Maxwell rheological model. For practical implementation in structural analyses, an approximate model for the isotropic homogenized fractured medium is formulated within the class of Burger models. Although the approximation is basically developed for short‐term and long‐term behaviors, numerical applications indicate that the approximate Burger model accurately reproduce the homogenized viscoelastic behavior also in the transient conditions.  相似文献   

Muon density measurements with the KASCADE central detector     

《Astroparticle Physics》2002,16(4):183-386

Frequency distributions of local muon densities in high-energy extensive air showers (EAS) are presented as signature of the primary cosmic ray energy spectrum in the knee region. Together with the gross shower variables like shower core position, angle of incidence, and the shower sizes, the KASCADE experiment is able to measure local muon densities for two different muon energy thresholds. The spectra have been reconstructed for various core distances, as well as for particular subsamples, classified on the basis of the shower size ratio N_μ/N_e. The measured density spectra of the total sample exhibit clear kinks reflecting the knee of the primary energy spectrum. While relatively sharp changes of the slopes are observed in the spectrum of EAS with small values of the shower size ratio, no such feature is detected at EAS of large N_μ/N_e ratio in the energy range of 1–10 PeV. Comparing the spectra for various thresholds and core distances with detailed Monte Carlo simulations the validity of EAS simulations is discussed.  相似文献   

关于高速滑坡的几个问题     

冯希杰  张骏 《地球科学与环境学报》1987,(4)

高速滑坡问题是当前工程地质和环境地质学界极为关注的问题。本文对诸如滑体高速滑动所对应的时空界限、滑体滑动过程所产生的热效应以及空气浮托力等高速滑坡产生机制问题进行了初步探讨。  相似文献